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Extinction of decomposable branching processes* 

Vatutin V.A.i Dyakonova E.E.^ 


Abstract 

The asymptotic behavior, as n —>■ oo of the conditional distribution 
of the number of particles in a decomposable critical branching process 
Z(m) = (Zi(m),^Ar(m)), with N types of particles at moment m = 
n — k, k = o(n), is investigated given that the extinction moment of the 
process is n. 
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1 Introduction 

We consider a Gallon-Watson branching process with N types of particles la¬ 
belled 1, 2,..., A in which a type i parent particle may produce children of types 
j > i only. Let rjij be the number of type j children produced by a type i parent 
particle. According to our assumption riij = 0 if z > j. 

In what follows we rather often use the offspring generating function of type 
N particles. For this reason, to simplify notation we put 

h{s) (1) 

Denote by the A-dimensional vector whose z-th component is equal to one 
while the remaining are equal to zero and let 0 = (0,..., 0) be an X-dimensional 
vector all whose components are equal to zero. 

Let 

Z(n) = (Zi(n),..., ZAr(rz)) 

be the population size at moment n G Z+ = {0,1,...} . We denote by 
mij{n) = E \Zj{n)\Z (0) = e^] 
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the expectations of the components of the vector Z(n). Let 

rriij = Wzj(l) = 

be the average number of the direct descendants of type j generated by a type i 
particle. 

We say that Hypothesis A is valid if the decomposable branching process 
with N types of particles is strongly critical, i.e. (see 0) 

= E = 1, i = (2) 

and, in addition, 

= E [r]i^^+l] e (0, oo), J = 1, 2,..., N -I, (3) 

E < 00, i = 1, k, j = i,i + 1,N, (4) 

with 

bi = ^Var[T]^^i] e (0,oo),i = 1,2,...,TV. (5) 

Thus, a particle of the process is able to produce the direct descendants of its own 
type, of the next in the order type, and (not necessarily, as direct descendants) 
of all the remaining in the order types, but not any preceding ones. 

In the sequel we assume (if otherwise is not stated) that Z(0) = ei, i.e. we 
suppose that the branching process in question is initiated at time n = 0 by a 
single particle of type 1. 

Denote by the extinction moment of the process. The aim of the present 
paper is to investigate the asymptotic behavior, as n —>■ oo of the conditional 
distribution of the number of particles in the decomposable critical branching 
process Z(to) = (Zi(to), ..., Zjv(ni)) with N types of particles at moment m = 
n — k,k = o(n) given Tjv = n. 

Decomposable branching processes of different structure have been investi¬ 
gated in a number of papers. 

We mention in this connection paper [5] dealing with the structure of the 
two-type decomposable critical Galton-Watson branching processes in which the 
total number of type 1 particles is fixed, and articles m-m in which, for A-type 
decomposable critical Markov branching processes asymptotic representations 
for the probability of the event {T^ > n} are found and Yaglom-type limit 
theorems are proved describing (under various restrictions) the distribution of 
the number of particles in these processes (and their reduced analogues) under 
the condition Tjv > n. 

In [3] the decomposable critical branching processes obeying the conditions 
of Hypothesis A were investigated under the assumption = n. In the men¬ 
tioned paper the conditional limit distributions are found for the properly scaled 
components of the vector 


Z(to) = {Zi{m),...,ZN{m)) 
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given that the parameter m = m(n) varies in such a way that 

TTl 

lim — = a; G [0,1). 

n—>-cxD 77, 

The following statement established in [3] is of a particular interest: 

Theorem 1 If Hypothesis A is valid and m = m(n) —>■ oo as n —>■ oo m such 
a way that m ^ xn, x G (0,1), then, for any Si G [0,1], i = 1, 2,iV — 1 and 
Xn > 0 


lim E 

n—¥C!0 


Zi(m) 

*1 


Zjv-i(m) 

■ ■ *Af-l 


exp 


— A AT 


ZN{m) 

bf^n 


Tn = n 


/ 1 + XMi-cr) ^ 

Vl + AAra;(l-a:)y (l + AArx(l 


a;)) 
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It follows from Theorem [T] that if Tjy = n and the parameter m = m(n) 
varies within the specified range, then the population consists (in the limit) of 
type N particles only. 

In the present paper, complimenting paper [3], we concentrate on the case 


hm — = 1. 

n—>-oo Ji 


Theorem 2 If Hypothesis A is valid and k = k(n) = n — m —>-00 as n ^ oo 
in such a way that k = o(n), then 


lim E 

n—¥oc 


^Zi(m) 


Ziv-i(m) 

^N-1 


exp < -Xn 


ZN{m) 

biqk 


Tn = n 


I 


(1 + Xn) 


Denote by hn{s) the n—th iteration of the probability generating function 
h{s). It is known (see, for instance, [1], page 93, formula (16)), that 


lim bNu'^ (hnis) - /i„(0)) = U{s) 

n—)-oo 


( 6 ) 


exists for any fixed s G [0,1), where U{s) is the generating function of the 
so-called harmonic measure. In addition, 


U{h{s)) = U{s) + l, sG[0,l). (7) 

The following theorem, complimenting Theorem [2] describes the final stage 
of the development of the process given its extinction at a distant moment n. 


Theorem 3 If Hypothesis A is valid and k = n — m = const as n ^ oo, then 


lim E 

n—¥C!0 


Zi{m) 

bi 


ZN-i{m) ZN{m) 

Sjv-l 


Tn = n 


Sn {U {sNhk+l{ 0 )) — U (sAr/lfe(0))) . 


Remark. We know by 0 that 


U {hk+i{0)) - U {hk{0)) = 1 


and, therefore, the limit distribution we have found in Theorem [3] is proper. 
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2 Auxiliary results 


We use the symbols Pi and Ei to denote the probability and expectation calcu¬ 
lated under the condition that a branching process is initiated at moment n = 0 
by a single particle of type i. Sometimes we write P and E for Pi and Ei, 
respectively. 

Introduce the constants 


Cn,N = 1/&W, 




1/2^ 


N-1 

n 




1/2^ 


i < N. (8) 


and 

Di = (6imi,i+i)^/^‘ci,i, * = 1,2, N. 
It is not difficult to check that 


Cl,AT = Djsf-i 



Dn-1 (cat,tv) 


1 / 2 ''-! 


Denote 


(9) 


( 10 ) 


Tki = min{n > 1 : Zk{n) + Zk+i{n) + ... -|- Zi{n) = 0|Z(0) = e/c} 

the extinction moment of the population consisting of particles of types k,k + 
1,..., *, given that the process was initiated at time r* = 0 by a single particle of 
type k. To simplify formulas we set = Tu. 

We fix > 2 and use, when it is needed, the notation 

70 = 0, 7,=7,(iV) = 2-("-0, i = l,2,...,N. 

The starting point of our arguments is the following theorem proved in [4] 
(see also [3): 

Theorem 4 Let Z(n), n = 0,1, ..., be a decomposable branching process meeting 
conditions m, m and 0- Then, as n ^ oo 

Pj(Z(n) 0) ^ Ci^Nn~^^'^ , (11) 

where Ci^N are the same as in 

This result was complemented in [3] by the following two statements the first 
of which is a local limit theorem. 

Theorem 5 (see m If Hypothesis A is valid, then, os r* —>■ oo 

i = l,2,...,N, (12) 

where 

9i,N = (13) 
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Corollary 6 (see W If Vn <^n, then 

lim P(Zi(0 + • • • + Zn-i{1) > 0\Tn = n)=0. (14) 

n—^oo 

Let rjrj {k, 1) be the number of type j daughter particles of the th particle 
of type r belonging to the fc—th generation and let 

N-l Tr Zr{k) 

= XI XI XI 

r=l k=0 q=l 

be the total number of type N daughter particles generated by all the particles 
of types 1, 2, TV — 1 ever born in the process given that the process is initiated 
at time n = 0 by a single particle of type 1. 

Asymptotic properties of the tail distribution of the random variable Wn 
are described in the next lemma. 


Lemma 7 (see Lemma 1). Let Hypothesis A be valid. Then, as 0 0 

l_E[e"®’^"|Z(0) = ei] (15) 

Proving the main results of the paper we will relay on the following state¬ 
ment. 


Lemma 8 If mN,N = 1, £ (0, oo) and k = k{n) = n — m ^ oo as n ^ oo 

in such a way that k = o{n), and 


s = exp 


A 

bNk 


,A > 0, 


then 

lim {hmis) - hmiO)) = 1. 

n—>-oo /c 

Proof. Let q = q{k) be a positive integer such that 


hq{0) <s< hq+l{0). 


Then 

hm+q{^) ~ h m (o)< ^ 771 ( 5 ) hm (o)< hm+q+l ( 0 ) /ijTT, (0). 

Clearly, q{k) —>■ 00 as fc —> 00 , and in view of the representation 

1 - hq{0){bNq)~^ , l-s^\{bNk)~^ 

and 

lim = 1, 

h^oo 1 — hqj^iyO) 
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we have q ^ kX~^, k —>■ oo. This and the local limit theorem for the critical 
Galton-Watson processes (see, for instance, [1], Corollary 1.9.1, p. 23) imply 
that, asl<CA:<Cm~n —>-00 


^ ^77^(0) - ^ ^ (hyrt-j-j-l-l (0) (0)) 

i=i 


E 


^ 6 Ar(m + j)^ bNn'^ hN^n^' 


Thus, 


bN^n^ 
lim sup — - - 

n—^oo 


{hm{s) - hmiO)) < 1. 


Similar arguments show that 


lim inf {h^{s) - /i„j(0)) > 1. 

n—^co /c 


The lemma is proved. 


Lemma 9 If Hypothesis A is valid, then, for any A > 0 


lim — 

n—J-oo fi^ 


7l 


■E 


Wn exp 


A— 

bisfU 


JlbNCl^N _ bNgi,N 
A^-T"! ~ A^-T'i 


Proof. Let 

WNik) = +wj^^ + ... + 

where the summands are independent random variables and = WnA = 
1,..., k. Setting 6 = \ (hNn)~^ in Lemma[7]and recalling (fT0)l . we obtain 


lim 

n—¥oo 


1 — exp 


A 

hMn 


Wn 



Cl,ArA'’'b 


Hence it follows that, for all A > 0 

A 


lim E 

n—^oo 


= lim E^' 

n—^oo 


exp 


bNn 


|Z(0) = K^]ei 


exp< - 


Inti 


W^MZ(0)=ei 


= expl-cpAfA"^^}. (16) 


Since the sequence of functions under the limit consists of analytical and uni¬ 
formly bounded functions in the domain {i?eA > 0} : 



A 

bNn 


Wn 


|Z(0) = ei 


< 1, 


it follows from the Montel theorem (see [6], Ch. VI, Section 7) that this sequence 
is compact. Moreover, since this sequence converges for real A > 0, the Vital! 
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theorem (see [B], Ch. VI, Section 8) and the uniqueness theorem for analytical 
functions imply convergence in (IT6l) for all A satisfying the condition {ReX > 0} . 
Moreover, according to the Weierstrass theorem (see [B], Ch. VI, Section 6) the 
derivatives of the prelimiting functions converge to the derivative of the limiting 
function in the specified domain. Whence, on account of the equality 


lim E 

n—>oo 


exp \ - f- |Z(0) = ei 


= 1 


it follows that 


1 


lim 

n^oo 


E 


Wat exp --— Wn 
onti 


El’ 


exp <{ - f- |Z(0) = ei 


= exp {-Cl,ATA"^^} 


or, in view of (ITBl) and (fTBl) 


lim 


jl-7l 


E 


Wat exp { --— Wn 
onti 


Ai-T'i ~ Ai-T"! 


The lemma is proved. 

Let 

Ik{m) = I{Zi{m) H-h Zkim) = 0} 

be the indicator of the event that there are no particles in of types 1, 2,..., k in 
the population at time m. We also agree to consider that Io{m) = 1. 


Corollary 10 If Hypothesis A is valid, then, for any A > 0 


lim 

n—>co 77 , 


l-7i 


E 


V^Arexp \ --^Wn [> lN-i{n^^^) 
' Onti 


bNgi,N 

Ai-T'i 


Proof. To check the validity of the statement of the lemma it is sufficient 
to note that, by virtue of (HU 

P (Zi(n2/3) + ... + Z^_i(n2/3) > o) = O = o . 

to make use of the equalities 

A 


E 

= E 

= E 


1 - exp -- Wn > ) lN-i{n'^^^) 


Inu 
A 


1 - exp --—ITw/Ar-i(n^/^) 


1 — exp < — 


bNn 

A 

bNn 


Wn 


- P 


(Zi(n2/3) + ... + Zjv-i(n2/3) >o) , 


and, replacing Wn by WNlN-iin"^^^), to repeat the arguments we have applied 
to prove Lemma m 
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Lemma 11 If Hypothesis A is valid, then, for any A > 0 and 1 <C fc <C n 


E 


Zjqfm) exp < —A 


r) l+7i 

n—¥cc) 


Proof. We consider the difference 


ZN{r 


A(m, fc; A) = E 
Clearly, 


exp < —A 


bj^k 


b^k 




bngi, 


N 


A2 


-E 


/{ZAr(m) =0}/Ar_i(n2/3) 


dA{m, k; A) 


1 


-E 


ZAr(m)exp |-A ^f 1 /Ar_i(n^/^) 


bNk 


dX bNk' 

Introduce, for m > Tn-i the quantity 

AT-l T,. Zr{k) 

n n n . 

r=l k=0 1=1 

Then, for < m < n 

Aim. t. A) = E [CH,„ (») - H,„ (0)) N-iin^'^'l 


bNk j 


where 

s = exp • 

Observe that, by the criticality condition 

hi+i(s) = h(hi(s)) > hi{s) 

and 

hi+i{s) - hi+i{0) = h{hi{s)) - h{hi{0)) < hi{s) - h/(0) 
for any fixed s G [0,1]. This yelds 

< ^ra (S) < . 

It is not difficult to check that for any s G [0,1] and < m < n 

^z„Mj^_i(^2/3)l _ E \l{ZN{m) = 0} J^_i(n2/3) 


E 

= E 


(H™ (s)-H™ (0))/Ar_i(n2/3) 


Using the inequalities 

tnA^b^U-. < J2r. 


1^1 


n 


1=1 

J 


bi 


i=i 


< 


< 


nd'-noi 


i=i 3=1 

[ai - bi) 


tn- 

Z=1 


(17) 


ai 


i=i 


1=1 



being valid for 0<b<bj<aj<a<l and nonnegative integers Vj, we obtain 


Wn (^m-TN-l(O))^" (^m(s) “ h^{0)) 

N-1 Tr Zr{k) / \ u /n\\ 

/j, ,\\b-m-k{s) — hm-kyo)) , . 

^ 

N-1 Tr. Zr{k) / \ u /n\\ 

Y^ Y^ /u ,\\'^m-k{s) — nm-k{'o)) , ^ 

^ - h — 7(7) -^ ^ 

r=l k=0 1=1 m ky j 

< fd^AT (/im(s))^'^ ^ {hm-TN-i{s) — hm-Tn-A^)) ■ 
Hence, on account of the condition Tjv-i < we conclude that 
(h^(s) - h^(0)) E \Wn (V-nV3(0))'^" lN-i{r7/^) 


< E 


{Hm ( S )- H ^( 0))/^_1 ( 71 ^/ 3 ) 


< (h^_„2/3(s) — /i^_„2/3(0)) E Wn {hm{s))^^ ^ lN-i{n^^^) 

Observe that if the parameter s has form dD, then for l<Cfc<Cm^n —>-00 

k 


hm(yS) hm(fi) hm—n^/si^s) /l^_„2/3(0) 


bN^V? 


in view of Lemma |8l and 


1 - hyn{s) - hm[0) 


bNn 


by the criticality condition. This, combined with Corollary 1101 yields 
A(m,fc;A) ^ {h„,{s) - h^m^\wN {hm{s))^^ lN-i{r7^^) 


bN^n? 

Thus, forl<CA:<Cm^n —>-00 


bNgi,Nn 


1-71 _ 


%i. 


N 




lim 


, 1+71 


n—foo k 

Let now A be a complex variable and 


A(m, fc; A) = 


gi.jv 

A ■ 


s = exp 


V 

bNk j 


ReX > 0. 


It is not difficult to check that for ReX > Aq > 0 and A: <C m ^ n —>■ 00 there 
exists a constant Ci(Ao) > 0 such that 


it 

l^m—n2/3('^) n2/3(0)| — t (^m—(I I) ^m—(0)) 


C-iCAo) 
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and, in addition, by Corollary (TU] there exists a constant C' 2 (Ao) > 0 such that 


t 7 i -1 


E 






Wn {hU\s\))'^^ < C2(Ao). 


Basing on these estimates and using the inequalities 

ri«?-rio?|=|i:w-!'r) ri d-rif? 

j=l j=l 1=1 j=l+l k=l 


1^1 


j=i^i 






being valid for 0 < \bj\ < |aj| < a < 1 and nonnegative integers Vj, we conclude 
that for ReX > Aq > 0 


,1+71 


^ |A(m,fc;A)| < - (/i„,_„v3(|s|) -/i,„_„2/a(0)) n^-^E /^_i( 

< CsiXo) 

for some constant C' 3 (Ao) > 0. This fact allows us, the same as in the proof of 
LemmaO to apply the Montel, Vital! and Weierstrass theorems and to conclude 
that for 1 <C fc <C n 


n—)-oo 


E 


-^N 


(m) exp I-A I 

d 6Arn^+'>'i 


dX k 


A(to, k; A) = 


bNgi,N 
A2 ■ 


The lemma is proved. 


3 Proofs of Theorems [2] and [3] 


Proof of Theorem By virtue of (HI it is sufficient to show that, for 
n ^ m ^ k ^ oo and k = n — m 


E 


exp -Ajv 1- /Ar-i(n^/^) 


b^k 


Tn = n 


1 


(Ajv + 1) 


2 • 


Put 


s = exp 


Ajv I 
bwk j ' 


Clearly, for m > 


E 


_^ZN(.m) 


lN-i(n^^^)] Tn = ' 


= E 




n2/3) 
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and, therefore, as /c —>■ oo 


s {hk+m - hfc(O)) E [Zjv(m) /Ar-i(n2/3) 


< E 




< s (hfc+i(0) - hfe(O)) E [zNim) (s/ife+i(0))^"(™)-' 

By the local limit theorem for the critical Galton-Watson processes 

1 


hk+iiO) - hfc(O) 


bNk'^ ’ 


oo. 


Further, 


shk+i (0) = exp I - + log hfe+i (0) | = exp | 


(1 + 0 ( 1 )) 


Using Lemma ITT] we conclude that for n ^ to ^ fc —>■ oo and k = n — m 

k'^bNgi.N 


E 


Zffim) exp |- (Aat + 1) (1 + o(l)) lN-i{n^^^) 


7j1+7i (A^ + 1) 


Since, as n —>■ oo 


in view of (HU, we have 


P {Tn = n) 


gi,N 

,^1+71 


E 




1 fc^fejvgi.Af 

6 Arfc 2 ^^2 

5l.Af „ , 1 

- 5 - ^ P (J AT = n)- TT 

7i1+7i (Aat + 1) (Aat + 1) 


as required. 

Proof of Theorem [3l The same as in the proof of Theoremj^l it is sufficient 
to show that, for a fixed k = n — m and n —>■ oo 


E 


^Z;v(m)j^_l(^2/3) 


Tn 



-)> s{U (s/ifc+i(0)) - U (s/ifc(0))). 


Put 

k* = [v^] j 'oi* = TO — [V^] = m — k* 

and write the representation 


E 

= E 


/ife. (s/i/c+i(0)))^"^™ ^ - (hfc. (shfc(O)))^"^" 
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Similarly to the estimates used in the proof of Theorem [5J we have 


s (hfe. (shk+iiO)) - hfc. (shfe(O))) E \zN{m*) (hfe. 


< E 




< s (hfc. {shk+im - hfe. (shfe(O))) E (hfe. (shfc+i(0)))^^(™*^-' 

According to (O 

lim {k*f (hfe. (shk+iiO)) - hk- (shfc(O))) = b~^ (C/(s/ifc+i(0)) - t7(shfc(0))). 

n—>oo 

Since 

1 - hk* (s/ife+i(0)) ^1- hk» (0) 


bi^k* 


as n —>■ 00 , and k* <C m* ^ n, then, according to Lemma ITT] 


lim 


,1+71 


n—)-oo ( h* 


{k* 


-E 


ZN{in*) exp < — 


ZNjrn*) 
bjsik* 




. = bjsfgi^N 


Thus, for any fixed k = n — m and n —>■ oo 

s iUishk+i{0)) - UishkiO))) bNgi,N{k*f 


E 


Tat = n 


bN{k 


*\2 


,1+71 


s {U{shk+iiO)) - U{shk{Q)))P {Tn = n), 


as required. 

Theorem [3] is proved. 
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